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Abstract 

Tropical rain forest has been a persistent feature in South America for at least 
55 million years. The future of the contemporary Amazon forest is uncertain, 
however, as the region is entering conditions with no past analogue, combining 
rapidly increasing air temperatures, high atmospheric carbon dioxide concentra- 
tions, possible extreme droughts, and extensive removal and modification by 
humans. Given the long-term Cenozoic cooling trend, it is unknown whether 
Amazon forests can tolerate air temperature increases, with suggestions that 
lowland forests lack warm-adapted taxa, leading to inevitable species losses. In 
response to this uncertainty, we posit a simple hypothesis: the older the age of 
a species prior to the Pleistocene, the warmer the climate it has previously sur- 
vived, with Pliocene (2.6-5 Ma) and late-Miocene (8-10 Ma) air temperature 
across Amazonia being similar to 2100 temperature projections under low and 
high carbon emission scenarios, respectively. Using comparative phylogeograph- 
ic analyses, we show that 9 of 12 widespread Amazon tree species have Pliocene 
or earlier lineages (>2.6 Ma), with seven dating from the Miocene (>5.6 Ma) 
and three >8 Ma. The remarkably old age of these species suggest that Amazon 
forests passed through warmth similar to 2100 levels and that, in the absence of 
other major environmental changes, near-term high temperature-induced mass 
species extinction is unlikely. 
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Introduction 

Warmth from anthropogenic greenhouse gas emissions 
may severely impact Amazon trees because the latitudinal 
distance species would need to move to maintain a con- 
stant temperature is large, implying that many species 
may be unable to track the future climate and become 
extinct (Colwell et al. 2008; Jones et al. 2009). More 
formally, climate envelope modeling arrives at a similar 
conclusion, but because upper-bound heat-tolerance in 
these models cannot exceed the observable maximum 
current surface air temperature, predicted extinctions 
may be inflated (Feeley and Silman 2010; Corlett 2011). 
Moreover, experiments show that tropical plants can 
photosynthesize and maintain a positive carbon balance 
under higher temperatures than those occurring today 



(Krause et al. 2010; Way and Oren 2010), thus suggesting 
that persistence may be possible for many lowland rain 
forest tree species. A complementary approach to assist 
in distinguishing the likely thermal tolerance of Amazon 
tree species is to consider historical evidence, because at 
times in the past, Amazon surface air temperatures have 
been higher than those today (Feeley and Silman 2010; 
Hoorn et al. 2010; Jaramillo et al. 2010; Haywood et al. 
2011). 

Some 56.3 Ma during the Paleocene-Eocene Thermal 
Maximum (PETM), global mean temperature increased 
by 5-6°C over a period of < 20 ka (Haywood et al. 
2011). Although many marine species became extinct, 
fossil pollen from the PETM showed an increase in 
tree diversity in three South American rainforest sites 
with abundant rainfall (Jaramillo et al. 2010). The mid- 
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Miocene climatic optimum (13-15 Ma) was the warmest 
time of the Neogene (the ca. 20 Ma Period encompassing 
the complete Pliocene and Miocene Epochs; (Zachos et al. 
2001). More recently, the early Pliocene (3.6-5.3 Ma) and 
parts of the Late Miocene (8-10 Ma) had surface air tem- 
peratures similar to IPCC predictions for the year 2100 
under low and high C0 2 emissions scenarios, respectively 
(IPCC 2007; You et al. 2009; Haywood et al. 2011). Fossil 
evidence suggests that megathermal (frost intolerant) for- 
ests were globally more extensive during the Neogene 
(Morley 2000) and possibly harbored higher levels of tree 
diversity than today's tropical forests (Hoorn et al. 2010). 
However, the relevance of Neogene warmth to the ther- 
mal tolerance of modern tropical trees depends, in part, 
on establishing the age of the tree species and their prior 
geographic distributions. If contemporary species extend 
back to the Pliocene or Late Miocene and did not have 
access to upper elevation thermal refugia, then they have 
likely survived air temperatures as high as or higher than 
those predicted for the coming decades. 

Here, we use a phylogeographic approach to determine 
if 12 contemporary Amazon tree species have endured 
historical air temperatures that exceed current thermal 
maxima. Because many of the 12 study species are rain 
forest specialists, their history provides proxy information 
about tropical rain forest distributions through time. 

Materials and Methods 

Study species and collection sites 

We obtained DNA sequences from the ITS and cpDNA 
genes and spacer regions from geographic samples of 12 
species (12 genera; 8 families) from the Amazon basin, 
Chocd region and Central America (Table 1; Fig. 1). The 
sampling sites were located on either side of the northern 
Andes Mountains, which have been a geographic barrier 
between Chocd and Central America (trans-Andes) and 
the Amazon basin (cis-Andes) through the Neogene 
(Hoorn et al. 2010), but more significantly so following 
the Merida uplift of the eastern Cordilleras at the Pliocene/ 
Pleistocene boundary, ca. 2.7 Ma (Helmens and van der 
Hammen 1994), which blocked a potential lowland dis- 
persal route along the Caribbean coast (Brumfield and 
Capparella 1996) (Fig. 1). 

Over 60% of the tree species (433 species) in Central 
Panama exhibit a cross-Andean disjunction and harbor 
Amazon basin populations (Dick et al. 2005). The 12 
study species were selected by comparing tree species lists 
from Center for Tropical Forest Sciences (CTFS) forest 
inventory plots in Barro Colorado Island, Panama (BCI) 
and Yasum' National Park, Ecuador (YAS). Trees were 
collected outside of the plots in both sites, and from 



remnant forests on the Pacific coast of Ecuador (Esmeral- 
das Province, ESM). Phylogeographic analyses were 
performed on >5 trees per species per site. Additional col- 
lections were obtained from herbarium specimens 
representing accessions from Costa Rica, Peru, Bolivia, 
Brazil, and French Guiana (Table SI). Voucher collections 
from ESM are stored in Quito (QCNE); Vouchers from 
YAS and BCI are stored at PUCE, MICH, and PAN 
herbaria. 

Species taxonomic boundaries have been delimited in 
recent monographs. The species have broad geographic 
ranges within the Neotropics, with a subset of species 
occurring in the West Indies and/or Africa. Andira iner- 
mis (W. Wright) H.B.K. (Fabaceae: Papilionoideae) is a 
canopy tree (height 10-15 m) of lowland forests in Africa 
and the Neotropics. Brosimum guianense (Aubl.) Huber 
(Moraceae) is a canopy tree (10-40 m) of mature terra 
firme and seasonally flooded forests. Ceiba pentandra (L.) 
Gaertn. (Malvaceae), kapok tree, is a wind-dispersed and 
bat-pollinated canopy emergent tree ( > 60 m) of dry and 
moist forests of the Neotropics and Africa (Dick et al. 
2007). Celtis schippii Standi. (Cannabaceae) is a canopy 
tree ( < 30 m) of mature lowland Neotropical forest. 
Chrysophyllum argenteum Jacq. (Sapotaceae) is a canopy 
tree (20-25 m) of lowland to premontane wet forests. 
Cupania cinerea Poepp. (Sapindaceae) is a mid-canopy 
tree (8-10 m) of upland and flooded forests in the conti- 
nental Neotropics and West Indies. Garcinia madruno 
(Kunth) Hammel (Clusiaceae) is a sub canopy tree 
(5-20 m) of mature moist forests. Ochroma pyramidale 
(Cav. ex Lam.) Urb. (Malvaceae), balsa, is a wind-dispersed 
pioneer of lowland and premontane forest. Palicourea 
guianensis Aubl. (Rubiaceae) is a clonal shrub or treelet 
(to about 3.5 m) of lowland and premontane moist forest 
and is often found at the edge of forest clearings. Poulse- 
nia armata (Miq.) Standi. (Moraceae) is a shade-tolerant 
subcanopy species found in relatively fertile soils and 
riparian habitats in moist lowland and premontane 
forests. Symphonia globulifera L. f. (Clusiaceae) is a shade- 
tolerant subcanopy to canopy tree ( > 30 meters) found 
along rivers and in upland forests in tropical Africa and 
the Neotropics including West Indies (Dick et al. 2003; 
Dick and Heuertz 2008). Tretna micrantha L. (Ulmaceae) 
is a small (2-20 m), bird-dispersed pioneer tree found 
along forest edges; it is a species complex that contains 
geographically overlapping clades (Yesson et al. 2004). 
This study is focused on a clade identified as Tretna sp. 
11 by Yesson et al. (2004). 

The study species are broadly representative of the 
Amazon tree flora: ten of the study species are animal- 
dispersed and two are wind-dispersed species, thus 
generally matching the proportions of seed dispersal 
syndromes found in tropical forests. Moreover, the 
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Table 1. Characteristics of study species, including geographic range (N = Neotropics, AA = Amphi-Atlantic (Africa + Neotropics), C = Caribbean 
Islands, At = Brazilian Atlantic Forests), successional stage (P = pioneer, GS = gap specialist; T = shade tolerant), and stature (C = canopy, 
U = understory, CE = canopy emergent) ('drought tolerant species), mean and standard error of tMRCA estimate, and 95% highest posterior 
density interval (HPD) (Bayesian analog of confidence intervals). 

Upper, Lower 95% 



Latin binomial 


Geographic Range 


Stage, Stature 


tMRCA (±SE) (Ma) 


HPD (Ma) 


Ceiba pentandra 


AA, C 


P, CE* 


0.2223 (0.00585) 


0.00014, 0.6975 


Andira inermis 


AA, C 


GS, C* 


1.4055 (0.00788) 


0.3713, 2.7062 


Trema micrantha 


N, C 


P, U* 


1.9167 (0.00905) 


0.53, 3.529 


Brosimum guianense 


N, At 


T, C 


2.9172 (0.01067) 


1.2462, 4.8233 


Palicourea guianensis 


N, C, At 


GS, U* 


3.0295 (0.01268) 


1.1007, 5.2625 


Cupania cinerea 


N 


T, C* 


6.4363 (0.01778 


3.3458, 9.5587 


Garcinia madruno 


N 


T, C 


6.9458 (0.01492) 


4.3041, 9.6291 


Celtis schippii 


N 


T C 


6.9848 (0.01663) 


4.0696, 10.0757 


Symphonia globulifera 


AA, C, At 


T, C 


7.0882 (0.01773) 


4.3201, 10.3423 


Ochroma pyramidale 


N, C 


P, C* 


8.0611 (0.01588) 


5.1609, 10.9799 


Poulsenia armata 


N 


T, C 


8.9379 (0.01882) 


5.58, 12.4139 


Chrysophyllum argenteum 


N, C 


T, C 


9.8766 (0.02049) 


6.4041, 13.9039 




Figure 1. Primary collection sites in (1) central Panama, (2) Western 
Ecuador (Esmeraldas Province), and (3) Amazonian Ecuador (Yasunf 
National Park). Additional collections were made in Brazil, Peru, 
French Guiana, and Bolivia for some species. The Andes and the 
llanos region presently form a strong geographic barrier between 
lowland moist forests, east and west of the Andes. The uplift of the 
Merida Cordillera occurred roughly at the Pliocene-Pleistocene 
boundary (ca. 2.7 Ma). 



species span the heliophile-shade tolerance spectrum and 
the understory, canopy and emergent stature continuum, 
two key axes of differentiation among lowland rain forest 
trees (Poorter et al. 2006) (Table 1). 

Laboratory methods 

Each DNA aliquot was assigned an identification number 
(Lab ID) linked to associated GPS coordinates, herbarium 
vouchers, or forest inventory tag numbers (Table SI). 
Molecular analyses were based on ca. 2400 bp of coding 



and non-coding chloroplast DNA (rfocL, trnH-psbA, psbB- 
psbV) and ca. 700 bp of the ITS region. Forward and 
reverse DNA strands were sequenced for each genomic 
region. Laboratory methods have described in detail 
elsewhere (Dick et al. 2007; Dick and Heuertz 2008). 
Haplotype networks were constructed using statistical 
parsimony (Clement et al. 2000). Homopolymer regions 
(cpDNA simple sequence repeats) were excluded from 
analyses due to DNA sequence alignment uncertainty and 
a high likelihood of homoplasy. 

Coalescent analysis 

We performed a Bayesian uncorrelated molecular clock 
analysis (Drummond et al. 2006) using the program 
BEAST (Drummond and Rambaut 2007) to estimate the 
time of most recent common ancestor (tMRCA) of the 
phylogeographic lineages for each species. The analysis 
was performed on the ITS dataset, for which many fossil- 
calibrated substitution rates are available (reviewed in 
Dick et al. (2003) and Kay et al. (2006)). Rate estimates 
for the ITS region in woody plants range from 0.5 to 
3.3 x 10~ 9 substitutions per site per year (s/s/yr) (eight 
studies) with a mean of 1.64 x 10~ 9 . We excluded rates 
for herbaceous plants, because of their twofold higher rate 
of nucleotide substitution (Smith and Donoghue 2008), 
and rates calibrated by vicariance events because, by 
excluding potential recent dispersal scenarios, these 
estimates are biased toward slow rates and older age 
estimates. We used the mean substitution rate, and the 
slowest and fastest rates for the lower and upper bounds 
of the standard deviation of the mean. The analyses were 
run using the HKY substitution model, a lognormal 
branch length distribution, and constant population size 
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for 10 7 generations, with data sampled every 10 
generations. Effective sample sizes (ESS) of coalescent 
genealogies exceeded 1000 for all tMCRA estimates. XML 
input files and MCMC log files are available upon 
request. 



Results 

Despite the slow mutation rates of the genetic markers, 
which are often used in phylogenetic studies, a majority 
of the species exhibited strong phylogeographic structure 
as indicated in the haplotype networks (Fig. 2). Mean 
tMRCA estimates ranged from 0.22 (±0.006) to 9.88 
(±0.021) Ma (Table 1). Nine of the 12 species had a 
Neogene tMRCA; seven species were of Miocene age 
(Table 1). Our results indicate that geographic lineages in 
most of the studied species (nine of 12 species) originated 
during the Neogene and that some of these populations 
passed through temperatures that were warmer than the 
present day. 

Discussion 

Although Pleistocene-age species have been documented 
within at least one large Amazon tree genus Inga (Rich- 
ardson et al. 2001), our study suggests that a cross-section 
of Amazon tree diversity may be much older. The results 
demonstrate the persistence of geographically structured 
populations of rain forest trees through several recent 
geological epochs and through major geological and cli- 
matic changes (Fig. 3). The nine Neogene-aged species 



have witnessed the Merida uplift, and the progressive 
closing of the Isthmus of Panama (ca. 3.1 Ma) (Coates 
and Obando 1996). The significant global climatic 
changes (Zachos et al. 2008; Willis et al. 2010) include 
(1) late-Miocene and early-Pliocene warm periods; (2) 
the intensification of Northern Hemisphere Glaciation 
(3.2-2.5 Ma); (3) the Mid-Pleistocene Revolution (ca. 
1-0.8 Ma), when mean tropical temperatures decreased, 
but variability increased on both orbital and millennial 
time-scales; and (4) the intensification of the Walker 
Circulation (ca. 1.9 Ma), which increased the seasonality 
of the hydrological cycle and the sensitivity of the tropics 
to precessional cycles (Maslin et al. 2005) (Fig. 3). Many 
of these older species occurred in regions in which no 
upper elevation refuges were available during periods of 
global warmth, suggesting that these species, and rain 
forest tree species more broadly, have higher thermal 
tolerance than niche models and some other studies have 
suggested (Thomas et al. 2004; Colwell et al. 2008). 
However, we cannot discount the possibly that some 
higher temperature tolerance may have been lost over time. 
This hypothesis should be a priority future research focus. 

Molecular clock analyses are subject to possible errors 
caused by rate variation across lineages, incorrect phylo- 
genetic placement of fossils for time calibration, and 
incomplete sorting of ancestral polymorphisms. Our study 
accounts for error associated with rate variation and 
genealogical stochasticity by incorporating the full range 
of relevant substitution rates into the Bayesian coalescent 
analysis. The application of even the fastest substitution 
rate in our study places the origin of 7 of the 12 species 
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Figure 2. Haplotype networks for nuclear ITS and chloroplast DNA sequence data. Circles represent haplotypes, colors are geographic locations; 
hatch marks represent mutational steps. ITS haplotypes are bounded by red rectangles. Additional geographic representation of some haplotypes 
is indicated by country abbreviations: Bolivia (Bo), Brazil (Br), French Guiana (FG), and Peru (Pe). 



) 2012 The Authors. Published by Blackwell Publishing Ltd. 



165 



Warmth Tolerance of Amazon Tree Species 



C.W. Dick et al. 



Benthic 
Warm 

3 



6180 (%») 

Cold^- 

4 5 



Global climate 



Regional events 



tMRCA of geographic populations (Ma) 




0.7-0.9 Ma 

Mid-Pleistocene 
Revolution 



2.0-1.8 Ma 

Intensification of 
Walker Circulation 



3.2-2.4 Ma 

Intensification of 
Northern Hemisphere 
Glaciation 



4.6-1 .8 Ma 

Progressive closure 
of Panama isthmus 



2.7 Ma 

Merida Andes 
uplift 



Glaciation of southern 
Greenland & Western 
Antarctic ice sheet 



I 



^ Terminal flood 

5.95-5.33 Ma 

Messinian 
Salinity Crisis 



Re-emergence and 
expansion of C4 plants 



Uplift of the 
Altiplano 



Zachos etal. (2008) 



0.0-0.7 • Ceiba 



0.4-2.7 # Andira 

0.5-3.5 • Trema 



1-1-5.3 Palicourea 

• Brosimum 

1.3-4.8 






3.4-9.6 # Cupania 

4.3-9.6 

4.3-10.3 
I4.1-irj.1» Symphonia 

Garcinia 

*" Celtis 

5.2-11.0 



5.6-12.4 Poulsenia 



6.4-13.9 %■ 



r .- Chrysophyllum 




Figure 3. tMRCA of geographic lineages and relevant corresponding surface air temperatures (Zachos et al. 2008) and climatic and geological 
events. Mean tMRCA estimates are shown in colored circles (Pleistocene blue, Pliocene green, Miocene red) accompanied by upper and lower 
95% HPD intervals. Large vertical arrows indicate period in which average surface air temperatures were warmer than present. Tree photographs 
were provided by Center for Tropical Forest Sciences. 



in the Neogene (Table 1). It is possible that some diver- 
gent haplotypes could be derived from introgression from 
related species. However, four of the older species are 
from genera (Symphonia, Celtis, Poulsenia, Ochroma) that 
each contains only one tree species in Neotropical forests, 
and thus these species are unable to hybridize or share 
ancestral variation with congeneric tree species. In fact, 
because we did not sample across the full geographic 
range of these species, some basal phylogeographic lin- 
eages may be unsampled. Our results, therefore, provide a 
conservative estimate of the species ages. 

The fossil record for our best-studied and most wide- 
spread species, Symphonia globulifera, corroborates the 
times scales presented here (Dick et al. 2003; Dick and 



Heuertz 2008). The earliest Symphonia pollen occurs in 
mid-Eocene deposits (ca. 45 Ma) in Nigeria, and S. globu- 
lifera pollen appears in mid-Miocene sediments in Vene- 
zuela and Brazil (ca. 15 Ma), and later in mid-Pliocene 
sediments in Mexico and Southeast Costa Rica. The mean 
tMRCA for S. globulifera reported here (ca. 7 Ma) is sub- 
stantially younger than its oldest Neotropical fossils (ca. 
15 Ma) (Dick et al. 2003). 

There are at least two demographic survival scenarios 
for these old species. Populations may have persisted in 
lowland forests suggesting that the lineages have survived 
past high air temperature levels. Alternatively, the species 
may have tracked higher elevation areas with tempera- 
tures matching their physiological optima. Populations 
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located in the western Amazon could find refuge on the 
slopes of the Andes, although these mountains were smal- 
ler in the Neogene (Hoorn et al. 2010). However, many 
of the species have disjunct populations in the West 
Indies and the Atlantic forests of Brazil (Table 1), where 
there were no nearby upper elevation refuges. Moreover, 
fossil evidence suggests forest stability in the core Amazon 
during the Neogene (Silva et al. 2010) rather than forest 
area contractions or the formation of lowland refugia. 

One possible difference between species survival 
through past warm climates and those projected for the 
21st century is the rate of temperature increase (Thomas 
et al. 2004; Botkin et al. 2007). Past climatic changes are 
thought to have included rapid warming events (Zachos 
et al. 2008) that, coupled with long tree generation times, 
suggest that acclimatization or standing adaptive variation 
was more important for climate response than major evo- 
lutionary innovation. Thus, decreases in effective popula- 
tion size and the loss of phylogeographic structure, such 
as that resulting from forest removal or fragmentation, 
are the causes for concern. The likely human-induced loss 
of genetic variation represents a novel difference when 
comparing 21st century Amazon rain forest with previous 
epochs, and may well impact the adaptive response of tree 
species to temperature increases and other environmental 
changes. 

The phylogeographic data presented here relate to spe- 
cies persistence and not to other forest attributes. 
Although experimental evidence suggests that tropical 
trees may physiologically tolerate a warmer climate than 
today (Krause et al. 2010), some experimental and eddy- 
covariance flux data demonstrate that tropical trees may 
be closed to, but have not exceeded, their temperature 
optima for carbon acquisition (Doughty and Goulden 
2008; Lloyd and Farquhar 2008; Lewis et al. 2009). All 
else being equal, higher temperatures in the future will 
increase carbon losses from respiration, while midday 
photosynthesis will be reduced, lowering productivity and 
tree biomass (Lloyd and Farquhar 2008; Lewis et al. 
2009). These impacts may then be reinforced by drought 
impacts (Lewis et al. 2011), or ameliorated by higher 
atmospheric carbon dioxide levels (Lloyd and Farquhar 
2008). Depending upon the strength of the fertilization 
effect of atmospheric carbon dioxide on tree growth, 
warming may reverse the current net carbon uptake by 
intact Amazon forests (Phillips et al. 2008), even if most 
species persist (e.g., see Galbraith et al. 2010; Rammig 
et al. 2010). 

Climate projections for Amazonia include not only 
higher air temperatures, but some models project an 
increasing frequency and intensity of drought events 
(Malhi et al. 2008), such as those seen in 2005 and 2010 
(Lewis et al. 2011). This may lead to lower biomass 



forest, replacement by novel forest types, or non-forest 
vegetation via forest dieback (Galbraith et al. 2010; Ram- 
mig et al. 2010). Yet, intriguingly, large-scale dieback dur- 
ing high temperature or drought events is not supported 
by paleoecological evidence (e.g., Mayle and Power 2008). 
The differences between model projections and paleoeco- 
logical studies may be due to the non-analogue condi- 
tions of 21st century Amazonia, but they could also stem 
from a lack of a paleoecological perspective within the 
biosphere modeling community. Our study cautiously 
suggests that many Amazon tree species can likely tolerate 
increasing air temperatures to some degree. This implies 
that drought events, direct human impacts (forest clear- 
ance, forest fragmentation, fires, and loss of seed dispers- 
ing animals), and their interactions may be more 
important immediate threats to the integrity of Amazon 
rain forests and therefore should remain a focus of con- 
servation policy. 
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online version of this article at the publisher's web site: 

Table SI. Accession and voucher information for 
sequenced individuals of the 12 study species. Samples 
were collected from localities within Ecuador (EC), Pan- 
ama (PA), Brazil (BR), Bolivia (Bo), and French Guiana 
(FG). Tag number references trees collected in permanent 
forest inventory plots. The ID refers to the catalog num- 
ber of DNA archived at the University of Michigan. 
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